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A theoretical study is performed on the entropy S a and the spin susceptibility Xs near the upper 
critical field H C 2 of s-wave type-II superconductors with arbitrary impurity concentrations. The 
changes of these quantities through H C 2 may be expressed as [S s (T, B)—S B (T, 0)] / [S n (T)—S s (T, 0)] = 
1 — c*s(l — B/H C 2) ~ (B I H C 2) as , for example, where B is the average flux density and S n denotes 
entropy in the normal state. It is found that the slopes as and a x at T = are identical, connected 
directly with the zero-energy density of states, and vary from 1.72 in the dirty limit to 0.5 ~ 0.6 in 
the clean limit. This mean-free-path dependence of as and a x at T = is quantitatively the same 
as that of the slope a p (T = 0) for the flux-flow resistivity studied previously. The result suggests 
that S S (B) and Xs(-B) near T — are convex downward (upward) in the dirty (clean) limit, deviating 
substantially from the linear behavior oc B/H C 2- The specific-heat jump at H C 2 also shows fairly 
large mean- free-path dependence. 

PACS numbers: 74.25.-q, 74.25. Op 



I. INTRODUCTION 

This paper considers the changes of the entropy S s and 
the spin susceptibility Xs through H C 2 for classic s-wave 
type-II superconductors. These quantities were calcu- 
lated by Makii^ in the dirty limit for superconducting 
alloys nearly 40 years ago. However, detailed studies on 
clean systems are still missing even for s-wave supercon- 
ductors. Writing these quantities as 
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the slopes as & n d a x will be obtained quantitatively for 
arbitrary impurity concentrations. The results near H C 2 
will also be useful for getting an insight into the behaviors 
over 0<B<H C 2- Indeed, a > 1 (a<l) indicates overall 
field dependence which is convex downward (upward), as 
realized from Eq. ([T|). 

It seems to have been widely accepted that various 
physical quantities of classic s-wave type-II superconduc- 
tors follow the linear field dependence with a = 1 at low 
temperatures. A theoretical basis for it is the density of 
states for a single vortex calculated by Caroli, de Gennes, 
and Matriconi2i4 However, few quantitative calculations 
have been carried out so far on the explicit field depen- 
dence. Recently, Ichioka et air performed a numerical 
study on the density of states of clean two-dimensional 
s-wave superconductors with k> 1 at T — 0.5T C . They 
found the exponent a = 0.67 for the overall field depen- 
dence of the zero-energy density of states. Also, experi- 
ments on the T-linear specific- heat coefficient 7 S (£>) for 
clean V 3 Si£ NbSe 2) 7 ' 8 ' 9 ' 10 and CeRu^ show marked up- 
ward deviations from the linear behavior "f n B/ H C 2- Even 
early experiments on "f s (B) for clean V and Nb indicate 



similar deviationsji2ii£ although not recognized explic- 
itly in those days due to the absence of a theory on clean 
systems. These results indicate that the field dependence 
with a < 1 is a general feature of clean s-wave supercon- 
ductors, as suggested by RamireziS 

Following the preceding works on the Maki para- 
meters^ and the flux-flow resistivity^ which will be re- 
ferred to as I and II, respectively, I here present a de- 
tailed study on S s and Xs near H C 2 at all temperatures. I 
thereby hope to clarify the k and mean- free-path (Z tr ) de- 
pendence of as and a x . Calculations are performed for 
both two and three dimensional isotropic systems to see 
the dependence of as and a x on detailed Fermi-surface 
structures. I also calculate the specific heat jump at H c2 
for various values of n and Ztr- To my knowledge, this 
kind of a systematic study has not been performed even 
for classic s-wave superconductors. 

Unlike the convention, I adopt the average flux density 
B in the bulk as an independent variable instead of the 
external field H . An advantage for it is that the irrele- 
vant region H < H c \ is automatically removed from the 
discussion on the field dependence. This distinction be- 
tween B and H becomes important for 1ow-k materials 
where H < H c \ occupies a substantial part of < H < H C 2 . 
Any experiment on the B dependence should be accom- 
panied by a careful measurement on the magnetization, 
especially for 1ow-k materials like Nb and V. 

Section [H] provides the formulation, Sec IIIII presents 
numerical results, and Sec. II VI summarizes the paper. I 
put /cb = 1 throughout. 



II. FORMULATION 

A. Entropy and Pauli paramagnetism 

As before^*^ I consider the s-wave pairing with an 
isotropic Fermi surface and s-wave impurity scattering 
in an external magnetic field H || z. The formulation 
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proceeds in exactly the same way for both the three di- 
mensional system and the two-dimensional system placed 
in the xy plane perpendicular to H. The vector potential 
in the bulk can be written ao 16 i 17 i 18 i 19 i 20 i 21 



A(r)=Bxy+ A(r) , 



(2) 



where B is the average flux density produced jointly by 
the external current and the supercurrent inside the sam- 
ple, and A expresses the spatially varying part of the 
magnetic field satisfying J V x A dr = 0. 

I first write down the expressions of the entropy and 
the magnetization in the presence of Pauli paramag- 
netism. As can be checked directly^ the effect can be 
included in the Eilenberger equations^ for the quasiclas- 
sical Green's functions /, and g by the replacement: 



e„ -> e' n = e„ - ifXB z- (V x A) , 



(3) 



where e n = (2n+l)TrT is the Matsubara energy and /i B 
is the Bohr magneton. The corresponding Eilenberger 
functional 2 ^ for the free-energy difference between the 
normal and superconducting states is given by 
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Here A is the pair potential, N(0) is the density of states 
per spin and per unit volume at the Fermi level, kp is 
the Fermi wavevector, and (• • • ) denotes Fermi-surface 
average satisfying (1) = 1. The quantity / is defined by^ 



I ee A*/ + A/ f + 2e' n [g- S gn(e n )} + h 



4r 



+h g(g) -i _ h P v F -df-fv F - d*p 



2t 



2[ ff +sgn(e„)] 



(5) 



where r is the relaxation time in the second-Born ap- 
proximation, vp is the Fermi velocity, and d denotes 
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we only have to differentiate with respect to the explicit 
temperature dependence in F. We thereby obtain 
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where S n = 2ir 2 N(0)VT/3 with V the volume of the sys- 
tem. In contrast, the expression of the external field H 
may be derived by applying the Doria-Gubernatis-Rainer 
scaling to Eq. (0J. 24 The details are given in Appendix A 
of I, and we obtain 

H = -4TrM nP + B + -L fdr (Vx A) 2 
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where M n p = 2fi%N(0)B denotes the normal-state mag- 
netization due to Pauli paramagnetism. We thus arrives 
at the expression of the magnetization from Pauli para- 
magnetism as 



MsP = MnP-i 
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When Pauli paramagnetism is negligible compared 
with the orbital diamagnetism by supercurrent, we can 
take the limit /x^^O in Eqs. J7J and © and retain only 
the leading-order terms. This results in e' n — > e n for Eq. 
1(7)1 . On the other hand, Eq. (0 is transformed by noting 
Eq. © into 



M sP = M nP 
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If the zero-field expression g — £ n /\/e 2 + |A| 2 is substi- 
tuted into Eq. iJTUjl with VxA = Bz, the terms in the 
square bracket reduces to the Yosida function^ 



The quasiclassical Green's functions / and g are con- 
nected by 3=(I-// t ) 1 / 2 sgn(e„) with P(e n , k F , r; /i B ) = 
f*(e n , — k F ,r; — // B ). The change of sign in ^ B is neces- 
sary here, because f = fulfil in the presence of Pauli 
paramagnetism. The functional derivatives of Eq. Q 
with respect to f\ A*, and A lead to the Eilenberger 
equation for /, the self-consistency equation for A(r), 
and the Maxwell equation for A, respectively. 

The expression of the entropy S s is obtained from Eq. 
by the thermodynamic relation: S s = S n — dF/dT. 
Considering the stationarity with respect to /, A, and A, 



B. Expressions near H C 2 

I now consider the cases where Pauli paramagnetism 
is small and provide explicit expressions to Eqs. (JJJ and 
HI near H C 2- From now on I adopt the units used 
previouslyi^i^ where the energy, the length, and the 
magnetic field are measured by the zero-temperature en- 
ergy gap A(0) at H = 0, the coherence length £ = 
?w F /A(0) with v F the Fermi velocity, and £?o = 0o/27r£ 2 
with (f>o ee hc/2e the flux quantum, respectively, with 
h=l. 
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First, /, g, and A are expanded up to the second order 
in A(r) as 



Differentiating Eq. Q 5|l with respect to T yields 



ff =(l-I/Wt / W) sgn(£ „) 
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Substituting them into Eqs. iJJJ and (|f U|) and using the 
Eilenberger equations for f 1 - 1 ' and .p 1 ^ to remove terms 
with vp-9, we obtain 



The quantity c?(/q ) /dH c2 has been calculated as Eqs. 
(31)-(32) of I to be 
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dr | A(r) | 2 — ^^(/^ 1 ^A A similar procedure enables us to obtain the expressions 
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Further, A(r) and /W near H c2 can be expanded in the 
basis functions V'Afq( r ) of the vortex lattice as 14 



A(r) = VFAoVo q (r), 



o!d{f^ ] )/de n and dffi/de n in Eq. flU as 



de n 



= -Y.(- l ) N (fNfN)^(s n ), (18a) 
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where is defined by Eq. (25) of I. Using Eqs. (JTSJ) 



(13a) and l|18a|) in Eq. (|14a|l . we obtain 



/«(£„, k F ,r) = VUA >T /«(£„, 6)e iN * ^ q (r) , 

(13b) 

where (0,(p) are the polar angles of vp with sin# — ► 1 
in two dimensions, iV denotes the Landau level, and q is 
an arbitrary chosen magnetic Bloch vector characterizing 
the broken translational symmetry of the flux- line lattice 
and specifying the core locations^ The coefficients Ao 
and fffl are both real for the relevant hexagonal lattice. 
Substituting these expressions into Eqs. I|12al) and l|12b|) 
and using the orthonormality of -0A'q( r ) and e lNv , we 
obtain 
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Except Aq(xH C 2—B, all the quantities in Eqs. (|14a|) and 
(|14b|l are to be evaluated at H c2 . 

It is possible to give an alternative expression to Eq. 
(|14afl using the equation for H c2 given by Eq. (33) of I: 
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Using Eq. H20|) we can also calculate the specific-heat 
jump at H c2 . It is given in conventional units as 1 
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where K2 is the Maki parameteni^ 2 ^ and /3a = 1.16. 
Equations JUJl and {2H> with Eqs. ifTTft. (JTBJ), and (J2U 

are the main analytic results of the paper. The quantities 
Ao, fjP, and K2 have been obtained in I. The explicit 
expression of fffi is given by 
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where K$ may be calculated efficiently by the procedure 
in Sec. IIF of I, with a change of definition of e n as 
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C. Analytic results at T = 



D. Analytic results in the dirty limit 



Now it will be shown that Eqs. I|14a|) and 114b|) reduce 
to an identical expression at T — for arbitrary impurity 
concentrations, which has the physical meaning of the 
zero-energy density of states. 

Let us start from Eq. I|14a|l where e n > and e n < 
yield the same contribution. Using this fact and Eq. 
(I18a.l l . it is transformed into 



I here summarize analytic results in the dirty limit 
r — > 0. First, the key quantities are calculated by 
choosing N cut = 1 in the procedure in Sec. IIF of I. The 
results are given by 
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The summation over n for T — > may be performed by 
using the Euler-Maclaurin formula and the asymptotic 
property (e„) ^e^ 1 (e n — >oo)m^ For example, 



Since (3 2 is of the order of 1/r, as shown below, 
may be approximated as (Kq) « (l/e n — f3 2 /i„} ~s n /(e 2 l + 
(/3 2 )). Using this (Kq) in Eq. and retaining only the 
leading-order contributions, we obtain 
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Equation 114bl) may be transformed similarly as 
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Thus, S s /S n = M s p/M n p, or equivalently, as 
T — for arbitrary impurity concentrations. 

Equations l|26a|) and (|26b|) have a simple physical 
meaning. Indeed, noting Eq. ifTTJl. I|13b|) . and Eq. (|18a|) . 

we find an alternative expression at T = 0: 
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which is nothing but the normalized density of states at 
e = 0. Thus, we have arrived at a simple result that 
the entropy and spin susceptibility at T = are both 
determined by the zero-energy density of states. 

The coefficient of Aq cx H c2 — B have been obtained 

in I. Also, / (1) '(0) in Eqs. (|2^ajl and l26bT) may be cal- 
culated efficiently from Eq. I|22|l by using the analytic 
expression^ 
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Notice that f± is smaller than /q by yfr. Substitution 
of Eq. i|3TI|) into Eq. lfT3|l leads to the equation for H c2 
obtained by Mak: 2 - and de Gennesi 2 ^ 

ln(T c /T)+V(l/2)- ip{x) = 0, (31) 

where ip is the digamma function, and x is defined by 

t((3 2 ) 1 , tH c2 
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with d = 2, 3 dimension of the system. As shown by 
Maki, 26 Eq. (|31fl can be solved near T = by using the 
asymptotic expression of ip(x) as 
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Thus f5 2 oc i? C 2 ~ r 1 , as assumed at the beginning. Dif- 
ferentiating Eq. (|31J) with respect to T, we obtain 
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Finally, k 2 and [An(-B)] 2 are calculated from Eqs. (34b) 
and (36) of I as 
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where kq is defined by Kn = <^o/27r£Q.ff c (0) with H c (0) 
the thermodynamic critical field at T = 0. Equation l|35|) 
agrees with the result by Caroli, Cyrot, and de Gennes. 28 
Now, let us substitute Eq. (fTTTTft into Eqs. (|14b(l and 
(|T§j) and use Eqs. (J2HJ). We thereby obtain 



with (3 = \/H c2 sin^/2v / 2. Hence, Eqs. 
at T = can be evaluated fairly easily. 
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In addition, Eq. I|18bjl is to be replaced by 
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Thus, M s p/M n p and S s /S n is the same at T = 0, in 
agreement with Eq. 1|27|) : they are both determined by 
the zero-energy density of states. Equation (|37b|l is the 
result obtained by MakL 2 Also, the expression 1 — 2A 2 , 
for the normalized zero-energy density of states at T = 
agrees with the result for the local density of states 
obtained by de Gennes4^i 

Equation ^ tells us that A 2 , = (1- B/H c2 )f3]^ as 
T — > for k 2 > 1. We hence find from Eqs. {p, (|57ajl . 
and (|37b|) that the initial slopes at T = for K2 ^ 1 are 
given by 



a s =a x = 2/ /3a = 1.72. 



(38) 



The results suggest the overall field dependence of S s and 
Xs at T — which is convex downward. Notice that the 
flux- flow resistivity pj at T = also has the same initial 
slope ol p — 1.72 in the dirty limit These results 

strongly suggest that the density of states at e = is 
mainly relevant to the physical properties of the vortex 
state at T=0. 
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Finally, the analytic results in the dirty limit are the 
same as those given in Sec. IID with a replacement of r 
by the transport life time Tt r defined through 
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E. The case with p-wave impurity scattering 



F. Numerical procedures 



If the p-wave impurity scattering is relevant, the fol- 
lowing additional terms appears on the right-hand side 
ofEq. ©: 

d /k.(k7t ) + ( ^).k / t +d g^ g )^ (3g) 



4r, 



2ti 



where (k'g) = (k'g(e„, k F , r)), for example, n is the p- 
wave relaxation time, and k is the unit vector along kp. 
However, Eqs. JZJ, ©, and ifTH)) remain unchanged once 
/ is modified as above. 

The corresponding calculations near H C 2 may be per- 
formed as described in Appendix A of I. It thereby follows 
that Eq. I|14|) and (|21|l are also valid together with Eqs. 

fTfl . DHJ, and iJUJl, where ffi is now given by 
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with 
L> = 
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I have adopted the same parameters as I and II to 
express different impurity concentrations: 



W^tr = l/27rT c r tr , l tr /l = r tr /r . 



(46) 



Numerical calculations of Eqs. I|14|l and (|21|l with Eqs. 
ljT7)l. ifTHj l. and ij2*U|) have been performed for each set of 
parameters by restricting every summation over the Mat- 
subara frequencies for those satisfying |e n | < e c . Choosing 
e c = 200 is sufficient to obtain an accuracy of ~0.1% for 
Eqs. I|14b|) and ||2TJ, whereas e c = 20000 (4000) is required 
for Eq. I|14al) in the dirty (clean) limit. Summations over 
Landau levels have been truncated at N = N cut where I 
put 7?-7V cut = 1 in the calculation of K$ > see Sec. IIF of 
I for the details. Enough convergence has been obtained 
by choosing A cut = 4, 40, 100, 200, 1500, and 4000 for 
£E/'tr = 50, 1.0, 0.5, 0.1, and 0.05, respectively. Finally, 
integrations over 9 have been performed by Simpson's 
formula with N cu t + 1 integration points for 0<9<ir/2. 



III. RESULTS 

Figures ^ an d [21 show temperature dependence of as 
and a x defined by Eqs. Ijlal) and ijlbj) . respectively, for 
different impurity concentrations parametrized by Eq. 
(|46p. They have been calculated in three dimensions for 
= 1.0 and «gl = 50. All the curves start from the 
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FIG. 1: Slope as as a function of T/T c for different impurity 
concentrations with d = 3, hr/l = l-0, and kgl = 50. 
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FIG. 3: Slope a(T = 0) = a s (T=0) = a x (T=0) as a function 
of ^E/hr for d = 2,3, h T /l = l,2, and kgl = 50. 
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FIG. 2: Slope a x as a function of T/T c for different impurity 
concentrations with d = 3, /tr/' = l-0, and kgl = 50. 



same value as = cv x = 0.862 at T = T C and develop differ- 
ences among different impurity concentrations at lower 
temperatures. The equality as — a x holds at T — 0, as 
shown by Eq. H27[) . and the value decreases from 1.72 
in the dirty limit to around 0.6 for ^/hi =0.1. Ac- 
cording to Eq. I|27|l , this variation in the slope at T = 
can be attributed to the dependence of the zero-energy 
density of states N s (0,B) upon the impurity concentra- 
tion. In particular, N s (0,B) in the dirty (clean) limit 
decreases more rapidly (mildly) than the linear behavior 
N n (0)B j H C 2 near H C 2- From this result, we expect the 
overall field dependence of the entropy and spin suscep- 
tibility at T — which is convex downward (upward) in 
the dirty (clean) limit, as realized from Eq. (JIJ. 

The difference between as and a x at finite tempera- 
tures is rather small, as expected from as = a x holding 
at T = and T c . In particular, the curves of as and a x 
in the dirty limit depend neither on the dimensions nor 
lti/l. However, the dependence develop gradually as the 
mean free path becomes longer. 

Figurc|21shows the slope a = as — a x at T = as a func- 
tion of ^E./lf[ for different combinations of dimensions 
and impurity scatterings. The four curves start from the 
same value 1.72 in the dirty limit, and decreases gradu- 



ally through unity towards 0.5 ~ 0.6 in the clean limit. 
However, we observe only small dependence of a(T = 0) 
on d and In /I- We thus realize that the zero-energy den- 
sity of states is mainly determined by the mean free path, 
and does not depend much on the dimensions nor the de- 
tails of the impurity scattering. 

It is interesting to note that the slope a p for the flux- 
flow resistivity p / , which was calculated previously^ for 
d — 2, 3, = 1.0, and kql = 50, show a complete nu- 
merical agreement at T = with the corresponding as 
and a x , i.e., as = a x = a p for arbitrary impurity concen- 
trations at T = 0. This fact indicates that pf at T = is 
also determined by the zero-energy density of states. 

Next, we examine the dependence of the slopes on the 
Ginzburg-Landau parameter kql ■ Figure 0] shows the 
same curves as Fig. ^ near the type-I-type-II boundary 
of kgl = 1. Each curve is shifted upwards from the corre- 
sponding one in Fig.^for kgl = 50, but the quantitative 
difference is rather small. This is also the case for a x . 
Thus, the slopes as and a x as a function of B do not 
have large kql dependence. 

Finally, Fig. [5] plots the specific-heat jump AC over 
T at H C 2 as a function of T/T c for different impurity 
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FIG. 4: Slope as as a function of T/T c for different impurity 
concentrations with d = 3, hi /I = 1.0, and «gl = 1. 



7 



1.0 r 




TIT 



FIG. 5: Specific-heat jump divided by T at H C 2, normalized 
by the corresponding quantity at T = T c and H = 0, as a 
function of T/T c with d = 3, Ztr/2 = 1, and kgl = 50. The 
curves are for ^E/^tr = oo, 10, 1, 0.25, and 0.1 from the top to 
the bottom. 



concentrations with d = 3, In /I = 1, and kgl = 50. It 
is normalized by the corresponding quantity at T = T c 
and F = 0, i.e., AC(T C )/T C = 1.43 in the weak-coupling 
model. The curves change gradually from almost T-lincar 
overall temperature dependence in the dirty limit to T 2 
dependence in the clean limit, and approach zero as ocT 2 
at lowest temperatures. 1 Although the ratio near T c is 
strongly dependent on kql asi£ 



AC(T cH )/T cH 

hm — : ; — r~, 

t cH ^t c AC(T C )/T C 



(2k 2 



GL 



1)& 



(47) 



the basic features pointed above are common among dif- 
ferent values of kgl, d=2,3, and ki/l— 1,2. 



all impurity concentrations. The results have been ex- 
pressed conveniently using the initial slopes as and a x 
defined by Eq. (Q. The main conclusions are summa- 
rized as follows: (i) as = a x holds both at T = and 
T = T c . (ii) as = a x = 0.862 at T = T C for all impurity 
concentrations, (iii) At T=0, the slope a decreases from 
1.72 in the dirty limit to 0.5^0.6 in the clean limit. This 
change is due completely to the mean-free-path depen- 
dence of the zero-energy density of states. The fact also 
suggests variation of the overall field dependence at T = 
from convex downward in the dirty limit to upward in the 
clean limit, (iv) The slopes have only small dependence 
on the dimensions and the details of the impurity scat- 
tering, (v) The slope a p of the flux-flow resistivity pf, 
which was calculated previously^ also shows a complete 
numerical agreement at T = with as and a x . This fact 
indicates that the zero-energy density of states is also 
responsible for pj at T = 0. 

The T-linear specific-heat coefficient J S (B) observed in 
clean materials£2i2ii2iiiii2ii£ show curves of a < 1 , which 
are in a qualitative agreement with the present calcu- 
lation. On the other hand, ~f s (B) for dirty samples 8,10 
follows well-accepted linear field dependence oc B/H C 2, 
which is apparently in contradiction with the present 
result in the dirty limit. It should be noted however 
that a careful experiment^ on pf shows field dependence 
near T = which is convex downward, and experimen- 
tally obtained a p agrees quantitatively with the dirty- 
limit theoryii£i22iSI Detailed experiments on the mean- 
free-path dependence of 7^ (B) and Pf(B) are desired to 
remove these discrepancies. 



IV. SUMMARY 

The entropy and the spin susceptibility near H C 2 have 
been calculated for s-wave type-II superconductors over 
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